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Abstract 

Interplanetary coronal mass ejections (ICMEs) have complex magnetic and density structures, which is the re- 
sult of their interaction with the structured solar wind and with previous eruptions. ICMEs are revealed by in situ 
measurements and in the past five years, through remote-sensing observations by heliospheric imagers. However, to 
understand and analyze these observations often requires the use of numerical modeling. It is because no instruments 
can yet provide a simple view of ICMEs in two or three dimensions. Numerical simulations can be used to determine 
the origin of a complex ejecta observed near Earth, or to analyze the origin, speed and extent of density structures 
observed remotely. Here, we review and discuss recent efforts to use numerical simulations of ICMEs to investigate 
the magnetic topology, density structure, energetics and kinematics of ICMEs in the interplanetary space. 

After reviewing existing numerical models of ICMEs, we first focus on numerical modeling in support of the SMEI 
and STEREO observations. 3-D simulations can help determining the origins of the fronts observed by SECCHI and 
SMEI, especially for complex events such as the January 24-25, 2007 CMEs. They can also be used to test different 
methods to derive ICME properties from remote observations, to predict and explain observational effects, and to 
understand the deceleration and deformation of ICMEs. In the last part, we focus on the numerical investigation of 
non-magnetic cloud ejecta. We discuss how simulations are crucial to understand the formation of non-twisted ejecta 
and the formation of complex ejecta due to the interaction of multiple ICMEs. 
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1. Introduction 

Coronal Mass Ejections (CMEs) have been observed 
re motely from space since th e early 1970s with OSO 



1995b : and (4) the COR-1 and COR-2 coronagraphs, 
part of the Sun-Earth Connection Co ronal and He 



liospheric I nvestigation (SECCHI, see \Howard et al. 



2002, 2008) onboard the two Solar-TErrestrial REla- 



7 {ToussA 1 19731) and Skylab jMacQueen et al.j \l223- Hons Observatory (STEREO) spacecraft {Kaiser et al 



Their association with geo-effective ejecta observed 
in-situ around Earth has be en made since the early 
1980s {Burlasa et al. , 1982 ). They are most often 
identified as moving bright fronts in coronagraphic 
instruments, such as: (1) the K-coronameters (cur- 
rently M k-4) of the Mauna Loa Solar Observatory 
(MLSO) {Elmore et a/.ll2003|): (2) t he co ronagraphs on- 
board the Solwind ( Sheelev et al.i 1980b and t he So - 



lar Maximum Mission (SMM) (MacOueen et al. , 1980); 



(3) the C2 and C3 coronagraphs, part of the Large 
Angle Solar COronagraph (LASCO) onboard the So- 
lar Heliospheric Observatory (SoHO) ( Brueckner et al.i 
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2008). These coronagraphs image the corona and 



the inner heliosphere up to heliospheric distances of 
about 0.1-0.15 AU (20-30 R Q ). Until the recent 
launch of Coriolis containing the Solar Mass Ejec- 
tion Imager (SMEI) and STEREO, there had been 
only a limited number of instruments dedicated to 
heliospheric observations of CMEs inside of 1 AU, 
mainly the in-situ and zodia cal polarimeters onboar d 
the two Helios spacecraft jjackson and Leinerk, 1985 ). 
To bridge the gap between the first 0.15 AU and 
Earth's orbit, numerical simulations of interplanetary 
coronal mass ejections (ICMEs) have been the pri- 
mary tool to investigate the propagation of ICMEs 
and their interaction with the interplanetary medium. 
The complexity of numerical models has been con- 
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Figure 1: Number density from different solar wind models (in cm -3 ). Upper left: Model of Roussev et al. (2003b) at 1 AU for the simulation of 
the August 24, 2002 CME, as reported m]Luggzj!lal] (2009a). Upper right: Model onCohen et all 1200711 at 1 AU for the simulation of the January 
24, 2007 CME, as reported m lLueaz. et a/.M2009cl) . Bottom: Model a ftOdstrcil etcd] 120051) at 0.1 AU for the May 12, 1997 CME. 



stantly increasing from simp le one-dimensional (1-D~ ) 
hydro - dynamic (HP) models (H undhau sen and Gentry. 



1969t \Gonzdlez-Esparza et all 120031) or m a gneto- 



hydrodynamics (MHD) models iDrver et all 11976 



Wu et all [2005) starting in the 1 970s to 2-D MHD 



models in the 1980s and 1990s ADrver et all 
Wuetal.i 19781: Vandasetal . 19951; \Rilev et al 



1979 



1997 



van der Hoist et all 120021) . More rece ntly, axi 



symmetric 3-D m o del ( s o-called 2.5-D) ( Wu et al\ 



1999; Carsill et al. , 2000; Schmidt and Careill, 2003) 



and fully 



Gombosi et al. 



120001; \Groth et all 120 00; Rile y et al. 
200 It Manchester et al. 2004b; Odstrcil et al.[ 2004 



3-D 

i r; 



MHD 



models 



Han et all 
I r 



1988 



Havashi 1 2005; Jaco bs et 312006; Naka mizo et al. 
20091 \Kataoka et all I2009F " have been developed and 
used to study ICMEs. Now, we can for the first time 
constrain and test these models with heliospheric re- 
mote observations of ICMEs. But, the use of simula- 
tions goes beyond reproducing observations. Numerical 
investigations of the propagation of ICMEs have been 
central to further our understanding of the deceleration 
of ICMEs in the heliosphere, the interaction of ICMEs 
with the bimodal solar wind and their associated defor- 
mation, the association of magnetic clouds observed in 
situ with CMEs and the interaction of multiple ICMEs 
to form complex ejecta. Most of these physical pro- 
cesses can not be investigated with the existing obser- 
vations and numerical simulations are required. In this 
review, we summarize past and present efforts to sim- 
ulate ICMEs in section 2, and we present some recent 



progresses regarding the simulation of ICMEs in sup- 
port of observations with the SMEI instrument and the 
STEREO mission in section 3. In section 4, we report 
on using simulations to test observational methods and 
we explore some of the major advancements in our un- 
derstanding of complex ejecta observed near Earth in 
section 5. The final conclusions are drawn in section 6. 

2. Numerical Simulations of ICMEs 

Heliospheric numerical models follow two main ap- 
proaches: (1) using an inner boundary past the critical 
point (where the fast Mach number is unity, usually be- 
yond 0.1 AU), or (2) coupling the heliospheric model 
self-consistently with a solar coronal model. Both ap- 
proaches have advantages and drawbacks. Mainly, be- 
cause the fast magnetosonic speed in the heliosphere 
is much less than that in the lower corona (usually by 
two orders of magnitude), the timestep in a heliospheric 
model is much greater than that in a solar coronal model 
(a fraction of a minute to a minute in an heliospheric 
model compared to about a tenth of a second in a coro- 
nal model). Therefore coupling an heliospheric and a 
solar coronal models (approach 2) is much more time 
intensive than just running an heliospheric model (ap- 
proach 1) and it can only be implemented on relatively 
large machines (64-256 CPUs). Because of the differ- 
ence in timesteps, a heliospheric model runs much faster 
than real time even on a small number of CPUs (4-32 
CPUs). A typical run requires only a few hundred of 
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CPU hours for a 0. 1 AU to 1 AU propagation. A cou- 
pled simulation can only run close to realtime on much 
larger machines (128-256 CPUs) and a typical run is 
5,000-50,000 CPU hours for a Sun to Earth propagation, 
depending on the resolution. However, determining the 
inner boundary conditions in approach 1 is a complex 
problem. It is especially true for the initiation of 3- 
D CMEs around 0. 1 AU with a realistic density struc- 
ture, magnetic structure and velocity field. By compari- 
son, initiating a CME on the solar surface (approach 2) 
yields complex and realistic structures at 0. 1 AU. More- 
over, approach 2 is the only one which can be used to 
study the interplanetary consequen ces of different ini- 
tiatio n mechanisms (see review by \Poedts and Jacobs , 
2010, in this volume). Here, we review some of the 
solar wind and CME models used for heliospheric nu- 
merical studies. 

2.1. Solar Wind Models 

The simplest numerical models of solar wind use an- 
alytical expressions for the solar wind speed, which in- 
creases approximatively like the logarithmic of the dis- 
tance, and for the solar wind density which decreases 
oc r~ 2 in the heliosphere. The next level of complex- 
ity is to define an axisymmetric bimodal solar wind, 
reflecting solar minimum conditions, where the slow 
solar wind has a speed of 300-400 km s and a den- 
sity of 5-15 crrT 3 at 1 AU and the fast wind a speed 
of 700-900 km s and a density of 1-5 crrT 3 . Such 
a solar wind can be set up using bot h approaches de 
scribed above, 
or i 



left an d right, show r esults f rom t he solar win d mod- 
els of \Roussev et al.\ d2003al) and \Cohen et al.\ (12007b 
based on a coupled corona-heliosphere code. The bot- 
to m panel shows the re sult from the heliospheric model 
o uOdstrcil et q/.l(l2005l) . All three models have a higher- 
density current sheet and lower density regions corre- 
sponding to the poles and coronal holes. The values 
obtained at 1 AU for the density between 3 and 20 cirT 3 
are typical to what is observed in-situ. 

2.2. Coronal Mass Ejection Models 

Most solar wind models are very similar in both ap- 
proaches, i.e. independently of the height where the in- 
ner boundary is set. In general, the same is not true 
for CME models. One exception is CME "models" 
which consist in adding a pressure or velocity pulse 
at the inner boundary. This simplest type of mod- 
els has been used for simulations starting close to the 
solar s urface \Usmanov and Dryer , 19951: Groth et al. 
12000b \Chane et al.l 120061) and pa st the fast magne 



tosonic point iOdstrcil et al. ,2002). CME models initi- 
ated at the solar surface hav e been r ecently reviewed in 
Roussev and Sokolov\ {2006) and in \Poedts and Jacobs] 



(2010). To summarize, there are two main classes of 
models for CMEs initiated at the solar surface: with 
or without a pre-existing flux rope. A pre-existing flux 
rope can be s uperposed out of equilibrium to the st eady - 

can 



rope can be s uperposed ou t or equilibrium to thestei 
state surface iLusaz. et 'all l2007t Wdth et al.l l2007h. 



It ha s been done in \Han et al.\ ([1988) 



emerge from below the ph otosphere (Amari et q/.[ |200 
Manchester et al. , 2004a : Archontis and Torokl |20o' 



Odstrcil and Pizz.o (1999) wi th an inner bounda ry at or can be superposed in equilibrium and erupt due to 



0.1 AU. and in Wu et al.\ i 19991). \Groth et al.\ d2000) and 
Gombosi et al] (2000) with an inner boundary at the so- 
lar surface. This class of models is particularly useful to 
study ICME deformation associated with the bimodal 
solar wind (see section 3.1). Finally, the most advanced 
class of models uses solar observations, typically mag- 
netogram, as input to create a full 3-D solar wind 
with realistic speeds and densities. Most models now 
use the Wang-Sheeley-Arge relation (W SA model, see 



flux cancellation o r boun d ary motions iWu et all 1999 



Chen and Shibata, 2000; Roussev et al.. 2003tl>. In 



the absence of a pre-exsiting flux rope, boundary mo- 
tions (shearing and/or converging motions) can result 



rope \Linker and Mikid 


19951: \Antiochos et al.l 


1999 


Amari et al.l 12000c \Linker et al.l l2003t \Lvnch et al. 


2005; 


van der Hoist et al.l 120071 \Aulanier et al. 


2010). 



Wans et al.l 1 1990c lAree and Pizza. 2000) between flux 



tube expansion and s peed at 1 AU. It can be used with a 
boundar y at 0. 1 AU iOdstrcil et al.l l2005h or at the solar 



surface iCohen et al.l 12007b \Oran et all 120091) . Other 
models such as thos e pr oposed by Usmanov] ( 1993b . 
Roussev et al\ d2003ab and \Odstrcil et al. ( 2004 ) involve 
empirical terms inspired by physical processes (wave- 
particle interaction) to create a distribution similar to 
that obtained from the WSA model. Figure Q] shows 
three examples of steady-state solar wind density dis- 
tributions used in recent simulations. The top panels, 



For heliospheric models with an inner boundary be- 
yond the critical point, these physical initiation mecha- 
nisms are very hard (for models using pre-existing out- 
of-equilibrium flux ropes) or even impossible (for mod- 
els invoking boundary motions) to adapt. An alterna- 
tive way, more realistic than simply using a pressure 
pulse, is to use geometrica l models of CMEs . It is often 



done with a cone model (Zhao et al. , 2002; Xieetal 



20041 12006b). which can be fitted to the particular char - 



acteristics of the studied event iOdstrcil et all [2004). 
Even using this model, however, it is hard to define 
the three-dimensional distribution of the density and 
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Figure 2: Upper left: Spheromak-type model for coronal mass ejecti ons ^Gibson and Low, 2000) as used in the m odel of\Manchester et al. I 2004b). 
Upper rig ht: Semi-circular flux rope model o ARoussev et al\ i2003al) used in the simulation o flLMgaz et al\ 420071) and based on the Titov-Demoulin 
flux rope (Titov and Demoulin, 1999). The image is zoomed on active region 9236 and shows the solar surface and an isosurface of d ensity color- 
coded with the radial mag netic field. 3-D magn etic field lines are drawn in white. Bottom: Coronal Mass Ejection (cone model, see Zh ao et q/.L 
\2Q0i\Odstrcil et a/.Ll2004l) used in the model o AOdstrcil and Pizz.d <2009h 



the magnetic field as well as the velocity field. Of- 
ten, the ICME is treated as a purely hydrodynami- 
cal perturbation (no magnetic field) and it is set with 
a uniform density and velocity inside the cone. Re- 
cently, it was improved to better reflect the flux rope 
shape of CMEs, t o use information obta ined from 
stereoscopic fitting ( Odstrcil and Pizza, 120091) and also 
to includ e non-uniform variat ion of the density and 
velocity (Hayashietal. Three examples of 



CME models used in recent heliospheric simulations 
are shown in Figure [2j illustrating two types of flux- 
rope models at the solar surface and one example of 
a cone model at 0.1 AU. The top panels show two 
flux ropes: a sp h eroma k-type flux rope developed by 
Gibson and Low (2000) (left) and a semi-cir cular flux 
rope developed by uitov and Demoulin ( 1 19991) (right), 
both added onto the solar surface. In both figures, the 
color shows the radial magnetic field strength in Gauss, 
illustrating how the flux ropes are added onto active 
regions. The b ottom panel shows the cone model of 



Zhao et al.\ d2002l) at the inner boundary of the code. 
The color shows the radial velocity, which illustrates 
how this ejection is added across the slow solar wind 



and extends into the fast solar wind from its onset. 

2.3. Reproducing Observations With Numerical Mod- 
els 

It is straightforward to extract a synthetic satellite 
file from a time-dependent numerical model, know- 
ing the spacecraft orbit. Such in situ measure- 
ments of the plasma parameters have been com 
pared to generic CMEs measur ed in situ hRilev et al. 



2003; Man chester et all l2004bl) and to specific events, 
recently the Halloween storm o f November 2003 

the November 24- 
20071). the M a y 1997 



(Toth etal, 2005; Wu et al. 



n '2005h 



25, 2000 CMEs (Luea z etal. 



CME iOdstrcil et all |2004 120051; \\Vu et all l2007b 



Cohen et al. , 20081) . Additionally, white-light images 



synthesis ha s been done i n heliospheric codes since the 



mid 2000s iLupaz et all l2005ah \Odstrcil et all 12005 
Sun et all 120081) . prior the launch of the STEREO mis 



sion. Some predictions of numerical models in sup- 
po rt of the STEREO missio n have been summarized 
in \Aschwanden et al. (120081) . The procedure to create 



these images is less straightforward than that to extract 
satellite data, but still relatively easy. Coronagraphs and 
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Figure 3: Top: CME interaction with the solar wind. Solid white lines streamline s display magnetic "streamlines" superimposed upon a color 
image of the veloci ty magnitude on the y — z plane (from [Manchester ~et al. [2004b)). Bottom: View of the November 14, 2007 CME reported by 
\Savani et al\i201$) at two different times in STEREO-B/COR-2 (left) and HI-1 (right ) fields-of-view. 



heliospheric imagers measure the Thomson-scattered 
light, which only depends on geometrical fa ctors and 
the electron density along a given line-of-sight &BillineA 
1966; Hundhausen, 1993). In a 3-D numerical simula- 
tion, the plasma density is known everywhere and a line 
integration can be used to calculate the Thomson inte- 
gral. Good agreement was found between observed and 
synthetic co r onagraphic line-of-sight i mages (e.g., see 
Linker et al. 1999; \Lusaz et all 1 2007) as well as with 



eclips e images (e.g., see Mik ic et ami9 99; Rusin et al. 



2010) . Some rece nt works (Ma nchester et all 2008 



Lusaz et all |2009c) have focused on making the syn- 



thetic procedure of creating coronagraphic images as 
close as possible from the real procedure, either by ac- 
counting fo r the F-corona (the light scattered by dust 
particles in \Manchester et al. , 2008) or by creating a 
background image for the MHD simulatio n using the 



proced ure developed for the real images iLusaz et al. 
20093). 



3. Support for Heliospheric Remote Imaging 

3.1. CME Deformation and Deceleration 

As the only way to study the 3-D properties of 
ICMEs, numerical simulations have been used to make 
predictions regarding the propagation of ICMEs in the 
interplanetary space. Numerical investigations have 
also helped deter mining which type of force s acts on 



ICMEs (e.g., see \Che 



ng which type oi torce s acts on 
A \\99&mippin\\2QQ& . ICMEs 



appear to be decelerated in the heliosphere by a drag- 
type force (see below) but there is also some evidence 
of a driving force, as some ICMEs a ccelerate or d o not 
decelerate much in the heliosphere. \Chen (1996) pro- 
posed that it is caused by the Lorentz force acting on 
the ICME. Alternatively, the presence of a fast stream 
ca tching up with an I CME is expected to accelerate 



it iGopalswamv et ail 120091) and it could be the rea- 
son why some CMEs do not decelerate. Often, obser- 
vations are not accurate enough to determine whether 
the ICMEs decelerate or not, and numerical simulations 
must then be used. 



In Lusaz et al.\ ( 2005 al) . using the simulation of 
Manchester et al. (2004b), we studied the evolution of 
the CME mass as the CME propagates in the interplane- 
tary space. We found that the mass increases by a factor 
of 5 from 1 R to 1 AU. This mass increase is associated 
with the solar wind material piled-up and compressed 
by the shock wave and incorporated into the CME. It is 
important to know this factor because, until the launch 
of STEREO and Coriolis, C ME masses could only be 
determined close to the Sun iVourlidas et all 120021) but 
mass appears to hav e some effects on th e associated 
geo-magnetic storms iFarrusia et al.luQQw . Numerical 
simulations are the easiest way to study how the CME 
mass changes as a CME propagates. 

A number of studies focused on the cause of the 
observed ICME heliospheric d e celeration (e . g., se e 
Riley et all 1 19971: \Carsili 12004 Wrsnak et al[ |2004 . 



Two physical models proposed to explain this deceler- 
ation are: (1) the "snow-plow" model where the mo- 
mentum of the system is conse rved and mas s increase 
is associated with deceleration uappini 120061) . and (2) 
the heliospheric drag acting on the CME to match its 
speed to that of the solar wind (without being associ- 
ated with mass increase). The drag force can be written 
C d (V C me - V sw )\Vcme - V sw \, where V C me and V sw are 
the CME and solar wind speeds, respectively. Estima- 
tions of the drag coefficient, C p, have been m ad e based 
on nu meri cal simulations in Vrsnak (2001), CarsilA 
d2004 and Wrsnak era/1(l2010h . among other works, re- 
porting values for the drag coefficient between 1 and 
10 AU~ l (except for extremely slow ICMEs), depend- 
ing on the initial speed and mass of the ICMEs, but re- 
maining relatively constant for a given CME from the 
solar corona to 1 AU. 

In all numerical studies of CMEs which include an 
idealized bimodal solar wind, it was found that wide 
ICMEs are significantly deformed as they propagate 
( Riley et al. , 1997 : Carsill and Schmidk 2002 ). It is due 
to the large drag force encountered near the current 
sheet, where the density is high and the solar wind speed 
is low. In contrast, ICMEs encounter a weaker drag 
force at high latitudes where the density is lower and the 
solar wind speed faster. The difference in the magnitude 
of the drag force can result in a concave-outward defor- 
mation of the ICME and its shock wave. Before the 
launch of the STEREO mission, testing findings from 
heliospheric numerical studies required statistical anal- 
yses u sing in-situ mea surements. Such an investiga- 



tion bw \Liu et al. (2008) found some signs that magnetic 
clouds during solar minimum conditions have concave- 
outward curvature, but it also found that the deflection 
flow ahead of the cloud was equatorward, which is the 



opposite behavior to that predicted by nume rical mod- 
els (for example, see [Manchester et q/.l 120051) . No sign 
of concav e-outward indenta t ion w as found in another 
study bv \Owens and Careili d2004 . While ICMEs ap- 
pear deformed and flattened in HI-1 field-of-view (see 
Figure 0, there has been only one clear instance of 
a change of indentation of the ICME leadi ng edge in 
SECCHI field-of-view iSavani et all l2010h . V-shape 



density en hancements are often o bserved at the back 
of CMEs iKahler and WebOi l2007h both by SMEI and 
SECCHI but this appears to be a different phenomenon. 
The deformation of a magnetic cloud and a ICME lead- 
ing edge requires a background solar wind which is 
at the same time steady and close to axi-symmetric. 
The presence of Corotating Interaction Regions (CIRs), 
dense streams and previous eruptions do not create con- 
ditions optimal to strongly deform ICMEs. The cur- 
rent solar minimum has h ad significantly mor e CIRs 
than the previous minimum AMason et al. , 2009), which 
might explain why such indentation changes have not 
been observed more often by SECCHI. It is also worth 
noting that recent simulations using non-axisymmetric 
solar wind models and/or multiple CMEs predict a flat- 
tening of the CME front but not nec essarily a change 



of its indentation (see simulations in Toth et ail 2007 



Lusaz et all 2007, 2009c; Web b et al. 



2009). 



3.2. Case Study of the January 24-25 CMEs 

After the launch of the two STEREO spacecraft in 
2006 October 26, one of the first CME events observed 
by the SECCHI suite was the two eruptions of January 
24 and 25, 2007 ^Harrison efa/.[l2008lE009l) . On 2007 
January 24 at 14:03UT, STEREO-A/CORl observed an 
eruption whose source region was behind the eastern 
limb of the Sun, with a speed of about 600-750 km s _1 . 
At 06:43UT on January 25, COR1 observed a sec- 
ond eruption also originating from behind the Sun and 
with a speed of about 1200-1400 km s . There was 
a ST£KEO-A/SECCHI data gap for 20 hours starting 
at 04:02UT on January 25, and a similar data gap for 
STEREO-B/SECCUl starting at 09:02UT on January 
25. Multiple fronts were detected in HI-2 after the data 
gap on January 26: the brightest one propagated from 
an elongation angle of 24.7° at 02:01UT to an angle of 
32.5° at 18:01UT (see Figure [4^). There are a num- 
ber of reasons to numerically study this particular series 
of events: (1) the 20-hour gap, which complicates the 
analysis of the data; (2) the fact that this is the first and 
still one of the best examples of CME-CME interaction 
observed by SECCHI; and, (3) the fact that the January 
25, 2007 was the fastest CME observed by SECCHI be- 
tween the launch of STEREO and 2010. 
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Figure 4: Comparison of synthetic (a) and real (b) J-maps for PA 69°, 
the center of the His image. The first tick on the x axis corresponds 
to 16:()0UT on January 24. Arrows point to similar tracks observed in 
the synthetic and real J-maps. 



3-D simu l ations o f these eruptions wer e reported in 
Lusaz et al.\ (E()08b). Lusazetal. (2009c), Web b et al. 



Lug 
(20( 



2009) and Odst rcil and Pizza (120091) . Here, we present 
a comparison between simulations and observations as 
well as important results from these simulations. Per- 
forming numerical investigations of specific events has 
several goals. First, it is the only true and appropriate 
test and validation for numerical models. Beyond sim- 
ply reproducing data, simulating real events is also one 
of the best ways to study the influence of the solar wind 
on ICMEs, by including not only an idealized bimodal 
solar wind, but a realistic 3-D solar wind, with dense 
streams, fast streams and CIRs. In the case of succes- 
sive CMEs, it is also important to move away from ide - 
alized simulations such as those of Lusaz et al. (2005b), 
Xions et al. d2006l) and \Xions et al.\ faOOll to integrate 
the complexity of real events and study the interaction 
between ICMEs of different speeds and directions and 
the structured solar wind. 



The simulation of \Lusaz et al.\ (I2009d) used the 



Space Weather Modeling Framework (SWMF, see 
Toth et al. , 2005b an d coupled a solar coronal model 



iCohen et all |2007) to an heliospheric model. The 



added onto the steady-state corona in the same active 
region behind the Sun's eastern limb, in a way similar 
to that shown in the top right panel of Figure [2] The 
two CMEs were initiated 16.67 hours apart and the 
second CME reaches the back of the first CME at 
approximatively 12:00UT on January 25. Later on, the 
overtaking shock wave exits the magnetic ejecta, and 
it then reaches the back of the dense sheath associated 
with the first CME around 18:00UT on January 25. The 
two shocks merge at around 06:00UT on January 26 
(see Figure 5 for 2-D cuts of the simulation at these 
different times). 

Figure [4] shows a comparison between synthetic and 
real time-elongation maps (J-maps) made from radial 
strips along the same position angle (PA): PA 69°, the 
apparent central position angle of the SECCHI instru- 
ments at the time of these ejections. The real J-map 
(Figure 4b) at this PA is complicated and hard to analyze 
due to the lack of observations on January 25. There 
are three tracks visible in HI-2 from the beginning of 
January 26 to about 14:00UT when a fourth track en- 
ters HI-2's field-of-view. From analytical models taking 
into consideration the initial speeds of the two ejections, 
it was expected that th ey would have me rged by the be- 
ginning of January 26 iWebb et al. [ l2009h . The presence 
of three tracks in the HI-2 field-of-view is therefore un- 
expected and hard to explain from the observations only. 
To fill the data gap and to understand the origin of the 
three tracks, we analyze the real J-map with the help of 
the J-map created from the simulation (Figure 4a). In 
the synthetic J-map at PA 69°, there are also three main 
distinct tracks (marked by black arrows). They corre- 
spond to: the first ejection, the second ejection and a 
dense stream from left to right, respectively. The dense 
stream is first compressed by the first CME and can be 
seen around 08:00UT on January 25 behind the first 
CME in Figure gji. At 04:00UT on January 26, there 
are three fronts past 20°, the brightest corresponding to 
the second CME. 

We believe that the three tracks in the real image cor- 
respond to the fronts identified in the synthetic image. 
According to our analysis, the third track corre s ponds t o 



the dense stream as explained in [Lusaz. et al. 



iponds to 
(2008b). 



CMEs were initiated w ith ou t- of-equilibrium fl ux ropes 
( Titov and Demouli 3 1 19991: Roussev et al. I H003b) 



It is visible before the interaction in the top left panel 
of Figure 5 at (X = 50 R Q , Y - 10 R e ) and after the 
interaction in the bottom middle panel of Figure 5 at 
(X = 20 R Q , Y = 60 R Q ). According to our 3-D simula- 
tion, this stream is present in the heliosphere prior to the 
CME passage, with a density about 50% larger than the 
average density at the same heliospheric distance. The 
simulation shows how the stream is compressed by the 
two CMEs, and due to the co-rotation, it enters the field- 



Figure 5: Left and Middle: View of the January 24 and 25, 2007 CMEs in the ecliptic plane (x,y plane) from the solar n orth (+z dir e ction) at four 
different times prior and during their interaction showing the density scaled by 1/R 2 as simulated with the SWMF in Lugaz et cd, 1 200') J). The 
white disk's radius is 16 R e and the black circle is the approximate projection of the Thomson sphere (relative to STEREO-A) onto the plane of 
the image. Earth's position is approximatively (-179, 1 10, -20)R e in this coordinate system. Right: View of the same CMEs from the solar north 
at two different times before and after the interaction as simulated with ENLIL in Webb et al. 1 2009). Earth's position is marked with the yellow 
circle at longitude 0° . 



of-view of HI-2 in January 26. The same dense stream 
is also clearly visible at b oth times in the simulation of 
Odstrcil and PizzAt00$) . 



The simulation of \Odstrcil and PizzA (120091) uses 
ENLIL and starts at 0.1 AU with boundary conditions 
given by the WSA model. The results are shown in the 
right column of Figure 5. For this simulation, the au- 
thors use a cone model input for the CMEs, which is 
purely hydrodynamic (no enhanced magnetic field in- 
side the ejecta). This simulation also predicts the merg- 
ing of the two ejections and the existence of a third 
dense front due to a corotating dense stream, but it 
predicts that this CIR will appear ahead of the CMEs 
in HI-2 field-of-view iOdstrcil and PizzA 120091) . The 
CMEs are found to merge around 06:00UT on January 
26 . One of the main differences with the simulation 
of Lusaz. et al. d2009cl) is the direct ion of propagation 



of the CMEs. In their simulations, \Odstrcil and Pizza 
(2009) initiate the CMEs on the eastern limb as seen 
from Earth. I mLueaz et a l. (2009c), the CMEs are initi- 
ated in active region 10940 which was about 20° behind 



the eastern limb at the time of the first eruption and 10° 
at the time of the second eruption. The start time of 
the ejection was also different between these two sim- 
ulations. In Odstrcil and PizzA d2009l) . the CMEs are 
initiated at 22 R Q at the time corre sponding to t he onse t 
time of the CMEs in LASCO. \n \Lusaz et al\ d2009ch . 
the time of the initiation of the CMEs was the time of 
first observation in COR1. 



The same events were reported in Web b et al. 
(2009) with a simulation with the Hakamada- 
Akas ofu-Fry Version 2 kinematic model (HAFV2 



see \Hakamada andAkasofuL Il982t \Frv et ail [2001). 
The steady-state solar wind is reconstructed using the 
Wang-Sheeley-Arge model with velocity pulses to ini- 
tiate the transients. The initial velocities are determined 
from the fit to the velocities derived from LASCO data 
and the duration of the pulse from the duration of the 
GOES X-ray profile. This model predicts that the two 
shocks have merged by 12:00UT on January 26. 



Thes e three heliospheric s imulat ions ALusaz et al 



(2009c|), \Odstrcil and PizzA d2009t) and \Webb etal 
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d2009h - are based on different numerical models and 
they use 3 different ways to initiate the CMEs: at the 
solar surface w ith an out-of-equilibr ium flux-rope (with 
the SMWF in \Lueaz et al.l l2009ch. at 2.5 wi th a 



velocity pulse (with HAFV2 in \Webb et al.l l2009h . at 



0. 1 AU with a non-magneti c cone model (with ENLIL 



in 



Odstrcil and Pizzo, 2009). They all predict that the 2 



CMEs are still distinct at 00:00U T on January 26 and 
hav e merged by 12:00UT. Bo th \Lueaz et all d2008bh 
and Odstrcil and PizzA d2009h show that one of the 
fronts observed in HI-2 is caused by a CI R or a dense 



stream , predicted to be behi nd the CMEs (Lusaz et al. 
2009ch or a head of them lOdstrcil and PizzA 12009). 



However, in \Lueaz etal. ( 2009d) . we also tried to ex 
plain the origin of another front observed in HI-2 at 
times when all three numerical models predict the 2 
CMEs have already merged (after 12:00UT on January 
26). We believe that the first bright front observed in HI- 
2 corresponds to the flanks of the first CME. The flanks 
of the two CMEs are distinct on January 26, although 
their noses have merged (see bottom left and bottom 
middle panels of Figure 5). Because the CMEs origi- 
nated from the limb, their flanks directly intersect the 
Thomson sphere and have a large contribution to the 
Thomson scattered signal. There are two main reasons 
why the noses and the flanks of the CMEs do not merge 
at the same time. First, the 2 CMEs propagate about 10° 
away from each other due to solar rotation. Second, the 
background solar wind is perturbed by the passage of 
the first CME, creating a more uniform background into 
which the second CME propagates. Consequently, the 
second CME flattens less than the first one (see bottom 
middle panel of Figure 5) and the distance between the 
CMEs' flanks is larger than that between the noses. 

Overall, the study of the January 24 and 25, 2007 
CMEs by three different groups show the importance 
of numerical studies in support of SECCHI observa- 
tions, especially for complex events, to correctly iden- 
tify the bright fronts observed and also, to learn more 
about CME-CME interaction. It is indeed impossible to 
understand the complex interaction of ICMEs without 
using a physics-based numerical code, especially in this 
case where there was a data gap in the observations. 

4. Testing Analyses Methods 

Magnetic clouds observed at Earth have been associ- 
ated with CMEs observed in the corona since the 1980s 



iBurlasa et al.l |1982|) . Associating complex structures 
observed in-situ with CMEs is not always as straight- 
forward. It is also not always clear from observations 



only how the 3-part structure of some CMEs (front, cav- 
ity, core) relates to the typical shock-sheath-magnetic 
cloud from in-situ measurem ents (see for example dis- 
Rilev et al.l 120081) . Numerical simulations 



cussion in 

of ICMEs are particularly important to study the ori- 
gin of the structures observed by in sit u instru ments on- 
board satellites. In Manchest er et al. d2004bl) . the au- 
thors compare the results at 1 AU of a 3-D MHD simu- 
lation initiated with a flux rope iGibson and Lom . 2000) 
with a typical three-part in-situ event (shock, sheath, 
magnetic cloud). They find good agreement between 
the simulation and the observation, which validates the 
results of the 3-D simulation. Moreover, synthetic line- 
of-sight images of the flux rope in the corona reveal that 
the flux rope is associated with the core and cavity of the 
CME, a result confirmed in a more recent study using 
a different initiation mechanism and solar wind model 



iRilev et all 120081) . This type of work has enabled sci- 
entists to associate the bright front of CMEs with piled- 
up mass and compressed material, corresponding ap- 
proximatively to the dense sheath observed in-situ prior 
to so me ICMEs. 

In utilev et al. (2004), the authors use a numerical 
simulation to blindly test different techniques to re- 
construct magnetic clouds from in-situ measurements. 
They initiate a CME at the solar surfac e via the "flux 
cancellation" method iRilev et al.l 12003) and follow the 
erupting flux rope all the way to 1 AU, where they 
extract a synthetic satellite measurement. Different 
magnetic cloud reconstruction methods are then tested. 
Numerical simulations can be used as a sort of con- 
trolled laboratory env ironment, in the following man- 
ner. For the study of utilev et al.\ (120044) . the magnetic 
topology of the magnetic cloud is exactly known, and 
the different reconstruction methods (force-free, Grad- 
Shafranov) can be rigorously tested, in particular with 
respect to the underlying assumptions. Theoretical er- 
ror analyses can be performed to quantify the accuracy 
of the reconstruction: for example, by choosing slightly 
different boundaries for the ejecta or a slightly different 
axes. However, it is nearly impossible to assess from 
observations only the errors associated with the assump- 
tions of the methods, for example the error associated 
with t he assumption of f orce-free or hydrostatic equilib- 
rium. \Rilev et al\ d2004h found that, when the spacecraft 
intersects the magnetic cloud off-axis, the resulting re- 
construction may not identify this ejecta as a magnetic 
cloud. This is an important result because only about 
30-40% of eject a observed at 1 AU are catalogued as 
magnetic clouds ^Burlasa et al. , 2003 ). From this study, 
it is now recognized that some of the remaining 60-70% 
are also magnetic clouds, but they are not recognized as 
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Time (hr) 



images with the position of the nose of the CME in 
the simulation. We also calculate the velocity from the 
reconstructed position and compare it with the CME 
speed from the simulation. This test of the approxima- 
tions is particularly interesting to do for the January 25, 
2007 CME because its speed varied significantly due 
to the interaction with the January 24 CME. We find 
(see Figure [6]l that the two approximations work well 
up to about 0.5 AU (here, corresponding to about 20-25 
hours) with errors in position of less than 10%. Beyond 
this distance, the Fixed-0 approximation over-estimates 
the CME position and predicts a CME acceleration (see 
right panel of Figure |6j. On contrary, the Point-P ap- 
proximation under-estimates the CME position and pre- 
dicts that the CME has decelerated down to the solar 
wind speed. These results do not reflect the 3-D sim- 
ulation (see right panel of Figure [6]) and the two ap- 
proximations should not be relied upon at large helio- 
spheric distances. Knowing precisely the position and 
kinematics of CME s is important to t est different drag 
models as done in Web b et al. d2009l) and to improve 
space weather forecasting. 



Figure 6: Position (top) and speed (bottom) of the January 25, 2007 
ICME at position angle 90° from the simulation along different radial 
trajectories, and, as derived from the synthetic SECCHI images with 
two different methods. 



such b ecause of the large impact parameter (Riley et al. 



2006). This is the case, for example, when the flank of a 
magnetic cloud passes over the spacecraft. In this case, 
some of the typical properties of magnetic clouds are 
not observed in the in-situ measurements. Other non- 
magnetic cloud ejecta are, however, intrinsically so and 
we discuss them in section 5. 

In \Lusaz et al. (2009b), we take a similar approach 
to test the two main techniques used to derive CME 
positions from elongation angle measurements made 
by SECCHI or SMEI. This wo rk is a continuation of 
the study by \Webb et al. \ (2009) and based on the nu- 
merical investigation in XLusaz et a l. (2009c). We ana- 
lyze the simulated data from synthetic SECCHI images 
with two approximations: th e Fixed-0 approximation 
&Kahler and VourlidaA 120051) . where the CME is as- 
sumed to be a single point propaga ting radially outward 
and th e Point-P approximation ( Vourlidas and Howarai 
2006), where the CME is assumed to be a spherical 
shell centered at the Sun. Because the analysis is based 
on a simulated CME, the position of the CME in the 
three-dimensional domain is known exactly. We com- 
pare the position derived from the synthetic line-of-sight 



5. Magnetic Clouds and Complex Ejecta 

A number of recent studies have focused on non- 
cloud ejecta, which are thought to be of three origins: 
(1) magnetic clouds (MCs) which are not recognized as 
such because they are observed off-axis (see previous 
section), (2) ejecta, which are originally MCs but that 
lose the characteristics of MCs due to interaction during 
their propagation, and, (3) ejecta, which are not initially 
MC-like. In this section, we review recent works on the 
last two types of ejecta. In-situ measurement at 1 AU 
provide all what is known about these ejecta. Therefore, 
numerical simulations are an invaluable tool to study the 
3-D magnetic topology of ejecta and the possible loss of 
the MC characteristics during propagation. 

5.1. CME-CME interaction 

With the improvement of coronagraphic observations 
and the presence of spacecraft measuring the solar wind 
in the outer heliosphere (Voyager, Ulysses), it was rec- 
ognized in the 1990s that successive ICMEs can merge 
with each other and also merge with CIRs to form a 
complex ejecta. In the outer heliosphere, these complex 
ejecta are sometimes referred to as a Merge d Interaction 



Regions (MIRs, see Burlasa et al. , 1997b . The same 



phenomenon also happens in the inner heliosphere, be- 
fore the CMEs reach Earth. One of the simplest cases 
is when two or more ejections interact but retain their 
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Figure 7: Left: ACE measurements of multiple-magnetic cloud event in March- April 2001 as reported for example in lVttoig et all (2 003^) The 
vertical lines show the position of the shock, the beginning and the end of the first cloud and the beginning and the end of the second cloud. Right: 
Results of our simulation of two successive CMEs as observed by a satellite at 1 AU as reported in Lugaz. et al. 1 2005b). The vertical lines show the 
position of the shock wave preceding the magnetic clouds, the beginning of the first cloud, a possible contact discontinuity inside the first cloud, 
the end of the first cloud and the beginning of the second cloud, from left to right respectively. Note the zone of enhanced /? between the two clouds 
in the observations and the simulation results. 



individual "identities" and form what is known as a 



multiple-magnetic cloud event iBurlasa et al.l 12002 



Wans et al.l 12002, 2003a). More complex cases hap- 
pen when it becomes nearly impossible to determine 
the solar origin of the complex ejecta observed in-situ. 
In this case, the characteristics of magnetic clouds are 



lost iBurlasa et al.l 2003; Wa ns et all 12004) . Numeri 



cal simulations are useful in both cases to determine the 
origin of the ejecta observed at Earth and the interaction 
processes which yield the observed complexity. It is 
possible to derive from coronagraphic observations the 
initial speed and delays between eruptions. From there, 
it is possible to infer (or guess) which eruptions might 
have had an impact on the form ation of a complex ejecta 
(e.g., see \Burlasa et a/.l l2003). However, it is more rig- 
orous to follow, via numerical simulations, successive 
CMEs from the solar surface to Earth's orbit, and to 
treat self-consistently their interaction. In the past five 



ye ars, such s i mulat ions have bee n performed in 1.5-D 
by \Wu et al.l d2005l) in 2.5-D by \Odstrcil et al\ d2005l) 



and 



Wu et al. (2002) a nd in 3- D bv\Schmidt and CarsilA 



d2004: \Lusazetal. (2005b) ; \Havashi et al.\ 72 006): 



Xion s et al. d2006il2007l):lLHgaz et all C2007 ; W u et all 
d2007ah : \Odstrcil and ftzzdd2009h . " 



Multiple-magnetic cloud events at 1 AU have the fol- 
lowing characteristics: (1) a uniform speed profile, (2) 
a shorter duration of the first cloud compare to the sec- 
ond one, and, (3) a r egion of high [3 between the clouds 
iWans et aZllixBa). The uniform speed appears to be 
associated with the ICME-ICME interaction, because 
the second ICME, in general, caught up with the first 
one and is expected to have a faster initial speed. It has 
been proposed that this uniform speed is due to the mo- 
mentum exchange associated with an elastic collision 
between the ICMEs, or with an inelasti c collision (as 
proposed by Farrusia and BerdichevsbX \2004) or with 
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Figure 8: Plasma speed in the y — z plane at four different times during the interaction of the two CMEs. "Streamlines" drawn in white illustrate 
the direction of the magnetic field in the plane. Top left: t = 13 h 45, as the trailing shock enters the first cloud. Top right: t = 19 h, as the trailing 
shock has just passed the center of the first cloud. Bottom left: t = 23 h 30, as the trailing shock enters the dense sheath of plasma associated with 
the leading shock. Bottom right t = 37 h. Note the uniformization of the speed in the two clouds at this time. 



the acceleration associated with the passage of a shock 
wave inside the first MC. These three mechanisms can 
also explain the shorter duration of the first cloud. With- 
out numerical simulation, it is i mpossible to rule o ut any 



of these three mechanisms. \a \Lueaz etal. (2005b), 



followed the interaction of two successive CMEs and 
we were able to reproduce and explain most of the fea- 
tures of multiple-magnetic cloud events (see Figure |7): 
the uniform speed profile, the high values of the pro- 
ton density and temperature in the sheath preceding the 
magnetic clouds and the existence of a region of high /3 
between the two clouds. We found that the speed of the 
event was relatively uniform, although the second CME 
caught up with the first one. This is due to a combina- 
tion of two factors. First, the shock driven by the faster, 
overtaking CME propagates through the first CME and 
accelerates it, it results in the speed of the two CMEs to 
be comparable. After the shock wave has already com- 
pressed the back of the first ejecta, the two ejecta collide 
which prevents the compressed first CME to expand and 
to slow down (see Figure[S}. We also confirmed that the 



region of high-/? between the two eruptions is associated 
with magnetic reconnection between the two clouds. 
The magnetic energy is transformed into thermal en- 
ergy, resulting in the enhanced /?. Last, we found that 
shock-shock merging may result in a very dense sheath 
preceding t he multiple-magnetic cloud event. It was 
proposed bv Warrueia et a l. (2006) that these extremely 
high densities can fill up the plasmaphere and it may be 
an essential pre-conditionment before the arrival of the 
large southward interplanetary magnetic field found in 
the ejecta to explain t he development of an intense ge- 
omagn etic st orm. In Wans et all ( 2003b ). Xione et al. 



(2006) and m \Lueaz et al.\ (12008a) . it was found that the 
passage of the overtaking shock wave inside the over- 
taken magnetic cloud indeed results in very large val- 
ues of the magnetic field. These simulations help ex- 
plain why instances of CME-CME interaction are asso- 
ciated with some of the lo ngest and largest geomagnetic 



storms jXie et a l., 2006a). 
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Figure 9: Left: Simulation results at t = 4 hrs showing a non-twisted ejecta resulting from shearing motions (model of Roussev et al. ( 2007)) applied 
to a dipolar active region. Right: Simulated satellite observations for this ejecta taken by an hypothetical spacecraft at 15 R G showing characteristics 
associated with a magnetic cloud and usually explained by a twisted flux rope (from Jacobs et al. (2002)). 



5.2. Non-Twisted Ejecta 

Other recent simulation efforts have focused on the 
existence of non-twisted ejecta. Magnetic clouds are 
commonly identified with twisted flux ropes. Some 
in-situ measurements of magnetic clouds appear to 
be in contradiction with this view. For example, in 
Most I et al. (2009b), the authors find that the amount of 



twist is maximal about 1/2 to 2/3rd of the way inside the 
ejecta and not on its outside, and the minimum is not on 
the ce ntral axis of the clo ud but near the outer bound- 
ary. In Mb stl et al. (2009a), the authors find a near con- 
stant twist in all the field lines of the magnetic cloud. 
The model of \Roussev et al\ (120071) attempts to mimic 
some of the shearing motions observed on the solar 
surface and relies on a realistic and complex magnetic 
topology using magnetogram data as input. This model 
was ap plied to the case stud ies of the April 21, 2002 
CMEs \Roussev et al. , 2007 ) as w ell as two successive 
CME s in August 22 and 24, 2002 hLueaz et aZ.Ll2009aL 



20101) . It was found that the CMEs do not develop into 
a twisted flux rope but, that, due to a series of recon- 
nection events, the ejected magnetic field lines exhibit 
writhe. Individual field lines are not twisted but have 
different orientations from one another, so that two con- 
secutiv e field lines have a slight angle between them- 
selves ( Tbrok et al. , 2010l) . The global aspect of the 
ejection exhibits smooth rotation of magnetic field lines 
without much twist (see Figure|5). The model was mod- 
ified and applied to simpler magnetic field topologies in 



Jacobs et al. (2009). The authors studied the shearing 
of a dipole embedded into a dipolar active region and 
into a quadrupolar active region. As the writhed mag- 
netic field lines propagate past an imaginary spacecraft, 
the satellite would observe a rotation of the magnetic 
field vector, reminiscent of a magnetic cloud, as can be 
seen in the bottom panel of Figure [9] Future work will 
test whethe r recons truction codes such as the one by 
Mbst l et al. (2009b) would mistake synthetic data pro- 
duced from simulations without twist for a twisted flux 
rope. If this is the case, the general paradigm of mag- 
netic cloud as a twisted flux rope would have to be re- 
vised. 

6. Discussions and Conclusion 

In this article, we have reviewed some of the recent 
simulations of ICME propagation and ICME-ICME in- 
teraction in the interplanetary space. We have discussed 
how simulations can be used in support of heliospheric 
missions (for example STEREO), and to study physi- 
cal questions such as shock-ICME interaction and helio- 
spheric drag. The two main types of heliopsheric MHD 
models have boundaries at the solar surface or beyond 
the critical point (typically 0.1 AU), and they can be 
used to study with different level of details ICME prop- 
agation and interaction. When these two different types 
of models, as well as a kinematic model, are used to 
simulate the same observed event (here, the January 24- 
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25, 2007 CMEs), they predict a similar behavior overall, 
but the differences in the details have important con- 
sequences. For example, the propagation of a shock 
wave inside a magnetic cloud has been found to greatly 
change the properties of the magn etic cloud as well a s 
the properties of the shock wave iLueaz. et al. . 2005b). 
Kinematics simulations and MHD simulations which do 
not use a magnetized ejecta are not expected to correctly 
reproduce the change in speed of a shock wave inside a 
magn etic cloud . This mig ht be the reason why the m od- 
els o flfrv et al. ( 200lb and lOdstrcil and Pizzo ( 20091) do 
not predict the CME-CME merging at the right helio- 
centric distance for this event. 

Numerical models have also been used to make pre- 
dictions of physical phenomena in the heliosphere. It 
is for example the case for the change of indenta- 
tion of ICMEs during solar minimum c onditions as 



they propagate close to the c urrent sheet ( Riley et al. 



2001; Manchester et al. , 2004b). While there have been 



some studies, which appeared to confirm these findings 



[Kah ler and "We bb. 2007: \Liu et al. . 2008; S avani et al. 



2010|), almost all CMEs observed by SECCHI do not 



show any sign of an indentation change. We suspect that 
the gre ater number o f CIR s during the current solar min- 



imum [ Mason et al. , 20091) compared to the previous 



one, is one potential cause for the discrepancy between 
simulations and observations. Simulations have also 
been used to study the nature of non-magnetic cloud 
ejecta, which are from 3 origins: (1) magnetic cloud- 
like ejecta which are not recogni zed as such be cause 



they are seen off-axis as shown m Utilev et al.\ (120041) . (2) 



complex eje cta resulting from th e interaction of succes- 



sive CMEs jLusaz. et al. , 2007b . and, (3) ejecta which 



are intrinsically non-twisted flux ropes, but which may 
have som e of the characteristi cs of magnetic clouds due 



to writhe iJacobs et a/.l l2009). 

Most heliospheric simulations of ICMEs use rel- 
atively simple models, from pressure or velocity 
pulses to out-of-equilibrium flux ropes initiated at 
the solar surface. By comparison, simulations of 
CMEs in the corona use more complex and re- 



emergence among other processes [Linker and Mikic 


1995; Antiochos et al.[ 


1999: \Amari et al.l 12003. 2004 


Manchester, 2007; Roussev et al. 


, 2007; Lynch et al. 


2008; 


van der Hoist et al., 


2009; Aulanier et al., 


2010) 



Most of these simulations follow the propagation of the 
CMEs up to few tens of solar radii at most and do not 
always have a realistic solar wind and/or coronal mag- 
netic field. With the continual improvement of com- 
putational resources and the existence of interplanetary 
measurements to compare simulations to, it is now im- 



portant to do heliospheric simulations of ICMEs using 
realistic models of CME initiation, to test these differ- 
ent models and study their consequences. Only in the 
past few years, th ey have been some attempts to do 
this, for example m \Lusaz et al. I d2010h with the Augus t 
24, 2002 CME and the model of \Roussev et al.\ d2007l) . 
On the other hand, simpler models ar e and will still be 



useful for space weather forecasting (LFrv et all 2003 



Taktakishvili et al. , 2009). In particular, they could be 



used to provide ensemble forecasting for space weather 
in a way comparable to what is done now with weather 
and climate modeling. Finally, combining heliospheric 
observations and the modeling of real events is ex- 
pected to improve our knowledge of CME-CME inter- 
action. Recent studies have explained why CME-CME 
interactions are associated with some of the largest ge- 
omagnetic storms: the compression of the magnetic 
field inside magnetic clouds and the formation of dense 
sheaths. It is important to study individual events to 
understand when and how these two phenomena hap- 
pen. We believe further modeling of real events will 
reveal the importance of studying individual events on 
case-by-case b asis n ot only at the Sun, as stated in 



Roussev et al.\ ([2007 ) but also in the heliosphere. We 
expect more progresses to be achieved in these three di- 
rections (testing of coronal models, ensemble forecast- 
ing of space weather and better understanding of indi- 
vidual events, in particular complex ones) during solar 
cycle 24. 
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